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We investigate an effective model Hamiltonian for organometaUic complexes that are widely used 
in optoelectronic devices. The two most important parameters in the model are J, the effective 
exchange interaction between the tt and tt* orbitals of the ligands, and e* , the renormalized energy 
gap between the highest occupied orbitals on the metal and on the ligand. We find that the degree of 
metal-to-ligand charge transfer (MLCT) character of the lowest triplet state is strongly dependent 
on the ratio e* / J. e* is purely a property of the complex and can be changed significantly by 
even small variations in the complex's chemistry, such as replacing substituents on the ligands. 
We find that that small changes in e* / J can cause large changes in the properties of the complex, 
including the lifetime of the triplet state and the probability of injected charges (electrons and holes) 
forming triplet excitations. These results give some insight into the observed large changes in the 
photophysical properties of organometaUic complexes caused by small changes in the ligands. 



I. INTRODUCTION 

OrganometaUic complexes are being developed as op- 
tically active materials for devices such as organic light 
emitting diodes (OLED),^"^ and organic photovoltaics 
(OP V) 1^"— The functionality of such complexes depends 
on their excited state properties. Of particular interest is 
the emission process, key to the function of organic light 
emitting diodes (OLED), 
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where C* (C) denotes the complex in its excited (ground) 
state. The reverse process, converting a photon into elec- 
tron and hole excitations, is the key process in OPV cells 
such as dye-sensitized solar cells (Gratzel cells) j^ To un- 
derstand and control these processes (in the design of new 
complexes) one needs to understand the relevant excited 
states of the complex. 

Two types of excited state typically dominate the be- 
havior of organometaUic materials (such as those illus- 
trated in Fig. [T|): hgand centered (LC) states and metal 
to ligand charge transfer (MLCT) states. Isolated lig- 
ands have a 'bright' high energy (UV) singlet transition 
(associated with a tt — >■ tt* excitation, equivalent to a 
LC transition). This transition has a molar absorptivity 
of order 10^ mol^^ cm~^fii. The corresponding tt — >■ tt* 
triplet transition has a much lower energy and is 'dark'.^^ 
When ligands are bonded to a transition metal to form a 
complex new singlet transitions in the visible region are 
observed (with molar absorptivity of order 10'^ mol~^ 
cm~^)jii These new features are typically attributed to 
MLCT transitions. The significant oscillator strength as- 
sociated with MLCT singlet states arises because of the 
hybridization of metal orbitals with tt* (unoccupied lig- 
and) orbitals, known as back-bonding j^^rJ^ 

The presence of the heavy transition metal ion pro- 
duces a spin-orbit interaction which mixes singlets and 
triplets with several important effectsJ^ First, this in- 
teraction allows states with dominant triplet character 
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FIG. 1. Examples of organometaUic complexes for opto- 
electronic applications. Both Pd(thpy)2 and Ir(ppy)3 are 
the basis of many complexes used in organic light emitting 
diodea^i^, while the benchmark N3 dya^ is used in organic 
dye-sensitized photovoltaic cells. 



to decay radiatively. Second, it also allows a rapid (^ 50 
fs) intersystem crossing (ISC) between singlet and triplet 
dominated statesJ^ In organic complexes with no transi- 
tion metal ion the spin-orbit coupling is much weaker, so 
triplet excitations tend to decay via non-radiative decay 
paths.iS-^ 

The character of the emitting state in organometaUic 
complexes has attracted considerable interest and debate. 
It is generally agreed to be predominantly triplet (due to 
its long radiative lifetime, which is associated with phos- 



phorescence) . However for a variety of complexes it has 
been labeled variously as LC, MLCT, or as a LC-MLCT 
hybrid.--^"^- For example, Li et al^ considered a family 
of OLED complexes and claimed that the emitting state 
is a superposition of a singlet MLCT state and a triplet 
LC state. Such spin hybridization is only possible due to 
spin-orbit coupling. In the case of very strong spin-orbit 
coupling, the predominantly triplet state can have a com- 
parable intensity to singlet emission^ Yersin et a/.^*'^^ 
claimed that the amount of zero-field splitting (ZFS) in 
the emitting triplet reflects the amount of MLCT char- 
acter in the state, and that there is "an empirical corre- 
lation between the amount of ZFS and the compound's 
potential for its use as emitter material in an OLED" ^ 
It has been suggested that the key effect of increasing 
MLCT character is to enhance the effect of spin-orbit 
coupling, which in turn increases the radiative rate of 
triplet excited states.— 

A key aspect to understanding the photophysical prop- 
erties of complex molecular materials is identifying the 
relevant frontier molecular orbitals and their interactions 
with one another. This allows one to define an effec- 
tive Hamiltonian which involves just a few parameters. 
Well known examples of this approach involve the Hiickel, 
Hubbard, Heisenberg and Pariser-Parr-Pople models for 
conjugated polyenes^^i With regard to organo-metallic 
complexes this approach has been applied to mixed va- 
lence binuclear systems including magnetic atoms in pro- 
teins (Hubbard and double exchange models)—, molecu- 
lar magnets,"^ and Anderson impurity models for cobalt 
based valence tautomers.'^^ 

Such semi-empirical approaches have significant ad- 
vantages which mean that they nicely complement first 
principles approaches such as DFT and high level ab ini- 
tio quantum chemistry.— First, effective Hamiltonians 
can help to reveal features that are common to a diverse 
range of materials. Second, since the number of degrees 
of freedom in the model is significantly less than in the 
actual material one does not necessarily have to make the 
approximations one would be forced to make if one did 
a complete quantum simulation of the actual material. 
For example, for large molecular systems one can also 
describe the nuclear dynamics quantum mechanically'^^ 
and include the effect of the environment such as the 
solvent i^i^ Such models can capture universal behav- 
ior that is not sensitive to microscopic details. For ex- 
ample, the single impurity Anderson model can describe 
the Kondo effect in a diverse range of systems including 
magnetic impurities in metals, quantum dots in semi- 
conductor heterostructures, carbon nanotubes, and sin- 
gle molecule transistors i^^*^ 

Here we develop and examine a model Hamiltonian 
which reproduces the key photo-physical properties of 
organometallic complexes. We apply the model to com- 
plexes in which the low energy ph ysic s is dominated by 
one ligand (see for example Refs. [l^ [2^, [SJ). We find 
that there are two key parameters in describing the low 
energy photophysics of these complexes - J, the spin- 



exchange interaction between ligand tt and tt* orbitals, 
and £*, the renormalized energy gap between the ligand 
HOMO and highest occupied metal orbital. We show 
that, through small changes to these parameters, single 
chemical substitutions can cause significant changes in 
the triplet excited state lifetime and the probability of 
triplet formation, causing large variations in efficiency. 
This paper is organized as follows. In Sec. II we in- 
troduce our model Hamiltonian, discussing its parame- 
ters and an appropriate basis in which to investigate its 
properties. In Sec. HI we analyze some approximate and 
exact solutions of the Hamiltonian, building an under- 
standing of the key features of the model. In Sec. IV we 
determine the effect of MLCT character on the radiative 
properties of the excited states. In Sec. V we study the 
probability of finding an excitation in the triplet state 
after charge injection. Sec. VI presents some concluding 
remarks. 



II. MODEL HAMILTONIAN 

In principle, one should determine an effective Hamil- 
tonian by integrating out high energy states. However, 
explicitly carrying out this procedure is prohibitively 
expensive for all but the smallest molecules i^ There- 
fore, in order to investigate correlation effects in these 
organometallic complexes we reduce the size of the ba- 
sis set in a way motivated by our understanding of the 
important physical processes. 

To this end we study the following model Hamiltonian: 

h = -jSh-Sl+ J2 {^^4Ma + alIl, 

- t« (HtM, + MlH„) - t^ [LlM, + mIl,) } 

+ UHriH-tnHi + ULnL-friLi + C/a/ nMt"-M; 

+ VHLnnnL + VnAinHnM + VMLnnnL, (1) 

where nn = Y^a^lHa, riL = Y.cr^l^'y^ "a/ = 

Y.aT.a'LlS<y<y'L'a, ^oa' = ((T^_„, , I7^_„, , <7^ ,,, ) whcrC 

(7^ ^, is the [a, a') element of the d Pauli spin matrix, a'^. 
In what follows it will be useful to have intuitive descrip- 
tions of the creation operators Hi, Ll^, and M^. Thus, 
we will refer to the state Hl\vac) (where \vac) is the state 
with (uh) = {um) = {til) = 0) as the 'HOMO level of 
the ligand', call ij^jwac) the 'LUMO level of the ligand', 
and call Ml\vac) the 'metal orbital'. It is important to 
stress, however, that the level Hl\vac) {Ll^\vac)) may 
well be very different from the highest occupied (lowest 
unoccupied) molecular orbital that one would find in a 
Hartree-Fock calculation. Similarly, Afjjwac) is likely to 
be very different from the atomic d orbital in the isolated 
transition metali^ As such one cannot describe the coeffi- 
cients of the Hamiltonian in terms of one electron orbitals 
like Uh = Jfdndr,rH^il)<PHt{l)y^rHi{^)^Hii2). 



Indeed, it is not possible, in general, to associate any 
static 'orbital distribution' with any of these states as or- 
bital relaxation will occur when the physical state of the 
complex changes. Rather, orbital relaxation (and higher 
order processes) are captured in the renornialization of 
the effective parameters of the modelj^*^ However, one 
expects that the model states correspond rather more 
closely to ones intuition of how a HOMO, LUMO, and 
metal orbital behave than the actual Hartree-Fock or- 
bitals do. Thus, this may provide new insights to the 
problems at hand»^ 

One might intuitively interpret the parameters as fol- 
lows: J is the effective exchange interaction between the 
ligand HOMO and the ligand LUMO, e is the effective 
difference in the energies of the metal orbital and the 
ligand HOMO, A is the effective HOMO-LUMO gap, 
t^ (t^) is the effective hopping amplitude between the 
metal orbital and the ligand HOMO (LUMO), Ui is the 
effective Coulomb repulsion between two electrons in or- 
bital i (where i € {i?, L, M}), and Vij is the effective 
Coulomb repulsion between an electron in orbital i and 
another in orbital j. This interpretation would seem 
to be equivalent to making an INDO (incomplete ne- 
glect of differential overlap) approximation^ after the 
physically motivated basis set reduction. However, we 
do not now make the Hartree-Fock approximation made 
in INDO calculations. Note that these parameters are 
renormalised from the values that would be found from 
direct computation. It is known that the parameter val- 
ues in small organic molecules, like the ligands consid- 
ered here (cf. Fig. [T]) can be significantly changed due to 
this renormalizatioui^Mirii Indeed, we will not attempt 
to calculate the values of the parameters here. Instead 
we construct a semi-empirical theory by comparing our 
model to experimental data in Appendix [X] 

In general, several ligands and metal orbitals could be 
involved in the low energy physics of the organometal- 
lic complexes we are investigating. For simplicity, here 
we focus on complexes in which one ligand dominates 
the low energy physics (for example Ru(NH3)2Cl2(bqdi) 
(Ref. 15), Ru(dcbpy)(bpy)2 (Ref. 37), or Pt(cnpmic) 
(Ref. ;26|)), interacting with only one fully occupied metal 
orbital, and leave the discussion of models with more 
metal and ligand sites for a later publication. In the non- 
interacting ground state {nn) = (nM) = 2 and (n^) = 0. 
Thus we consider a three site model with four electrons. 



Basis States 

One expects that the LUMO will be much higher in 
energy than the metal state (A — e :^ t^, see Appendix). 
Therefore, the effect of t^ on the eigenstates will be small 
compared to the other interactions. As such n^ is almost 
a quantum number so it is convenient to work in a basis 
of states which are eigenstates of ul- 

We define a 'reference state' for this model 



the TT-L = state, with energy Eq = {0\H\0) = 2e + UH + 
Um + ^Vhm, and where \vac) is the vacuum state. 

We will also choose our basis states to be eigenstates of 
the exchange interaction, J, that is, singlets and triplets. 
We define the til = 1 metal to ligand charge transfer 
(MLCT) triplet states as 

\HlLCT') = {4a/J0), 4^/^10), 
1 



^A,f._ 4^4)10)}, 



and the ligand centered (LC) triplet states as 

\'Lc') ^ |4i/^io), 4fftio), -^(4i/t-4^i)io)}, 

where the prefix "^ indicates the spin degeneracy and the 
suffix ^ is the value of n^ in these states (since these are 
triplets between an electron and a hole in the reference 
state, the S^ = states appear to have the opposite 
phase relation to usual). Similarly, the MLCT singlet is 

\HlLCT') = ^{l\m^ + lIm^)\0), 

and the LC singlet is 

\'lc')^^{l\h^ + lIh^)\o). 

When we go to the n^ = 2 states we have three singlets, 
the MLCT 



theLC 



l^MLCT^) = l\lIMiM^\0), 



\^lc^) = l\lIhiH^\o), 



and the metal-HOMO (MH) singlet 

and only one kind of triplet, 

\H'IH^) = {44i/^A4|0), L\LlH^AmO), 



-^l\lI{H^M^ + H^M^m}. 



|0) 



iJ.^ff|M|Af||uac) 
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When examining the results of our model, it will be 
helpful to consider numerical results for a particular 
set of parameter values. In the Appendix we estimate 
parameter values relevant to complexes of interest for 
optoelectronic applications (such as those in Fig. [1}. 
We use the following parameter set as a typical exam- 
ple of these values: J = 1 eV, A = 3 eV, i^ = 0.1 
eV, t^ = 0.1 eV, Um ^ Uh ^ Ul ^ U ^ 3 eV, 
Vhl = Vhm = Vlm = V = i gM (note that the approx- 
imation here that Uh — Vhl is physically reasonable as 
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Eigenstates 

FIG. 2. The eigenstates of the singlet part of the model 
Hamiltonian, Eq. [T] from lowest energy (1) to highest (6), in 
terms of singlet basis states which are eigenstates of ul , using 
the example parameter set J = 1 eV, A = 3 eV, e — 0.25 eV, 
t" = 0.1 eV, t^ = 0.1 eV, Um = Uh = Ul = U ^ 3 eV, 
Vhl ~ Vhm = Vlm = V = 3 eV (see Appendix for more de- 
tails). It is clear that for this set of parameters the eigenstates 
are predominantly one basis state, and that the coupling be- 
tween blocks of different hl is a small effect. These results 
hold for physically reasonable parameter values. See Supp. 
Info, for the equivalent figure at different parameter values.— 
Figs. S1-S6 in the Supp. Info, show that this picture does not 
break down except near (avoided) level crossings, and that the 
first excited state has a large overlap with the MLCT^ basis 
state for a wide range of parameters around our typical set. 



the H^vac) and L^vac) states are delocalised and have 
a large spatial overlap; a detailed derivation of this result 
is given in the Appendix). It is not possible to estimate 
e (the energy gap between the HOMO and the metal 
orbital) from what is known experimentally about the 
ligand as it is, intrinsically, a property of the complex. 
Indeed, we will show below that many of the important 
properties of the complex depend sensitively on s. For 
consistency and concreteness we take £ = 0.5 eV when 
we wish to investigate how the properties of the complex 
vary with another parameter. In the Supplementary In- 
formation we show that our main conclusions are robust 
to variations of these parameters in physically reasonable 
ranges and present calculations for many alternative pa- 
rameter sets. 

Fig. [5]shows that for this typical set of parameters that 
ul is indeed almost a quantum number for the singlet 
states, particularly in the low energy singlets. 

It is useful to write the Hamiltonian matrix in the basis 
defined above. We can write the singlet sector of the 
Hamiltonian in terms of blocks of constant n^, H^'^ , and 

the matrices that couple them Tg ^ ' ^ , 
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and similarly for the triplet sector of the Hamiltonian, 
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where T^"'""*=T;'"'"". 



The only ul = basis state is the reference state, H'^ — (2e + Uh + Um + 4Vhm). The singlet iil = 1 block is 



Hl = 



A + e + UH + 2Vhl + 2Vhm + Vlm 






f-ff 



t" A + 2e + ii + Um + Vhl + 2Vhm + 2V^m , 

where the states above the matrix indicate the basis in which we write the matrix. The triplet n^ = 1 block is 



(4) 



Hrp 



l^MLCT^) 
A + e + UH + 2Vhl + 2Vhm + Vlm 



H 



\^LC^) 



H 



t 

A + 2e -i + UM + VHL + 2Vhm + 2Vlm 



(5) 



the singlet ul — 2 block is 



Hi 



2A + Uh + Ul + Whl V2t" 

V2t" 2A + e + UL+2VHL + VHM + 2VLM V2t" \, (6) 

V2t" 2A + 2e + UL + Um + Wlm 



r 



and there is only one state in the triplet n^ = 2 block, the The coupling matrices are 
\H'IH^) state, H^ = {2A+e+UL+2VHL+VHM+2VLM). 
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amongst the singlet states, and for the triplets 
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III. 



APPROXIMATE AND EXACT SOLUTIONS 
OF THE MODEL HAMILTONIAN 



We will now investigate both the analytical solutions 
of Eq. [1] in the limit of t^ = 0,t" ^ and the exact 
solutions numerically, i.e., we find the full configuration 
interaction solution, to gain some insight into the behav- 
ior of the model. 

A typical value of A is around 3 eV, much larger than 
the typical i^ values of 0.1 eV (see Appendix). This 
means that the hybridization between states of different 
ul will be small compared to the hybridization between 
states of the same til, as is clear from Fig. [51 As such, we 
expect the approximation that i''^ = is close to the exact 
solution. In the limit i^ = all of the T"^'"'^ matrices 
are null, and hence ul is a good quantum number. 

The lowest excited states will come from the til — 1 
subspace, provided A is large, as we expect it to be in 
complexes with OLED and OPV applications (see the 
Appendix) . The eigenvalues of this subspace are 
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(i)± 
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[3-4§]^ 



+ 6Vhl + 2Vhm 



2Um + ^Uh 



(£* + [3-4§]-) +4(t^)2 
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where 



e*=e + UM-UH + Vl 



M 



Vl 



HL 



(7) 



(8) 



is the effective HOMO-metal energy gap relevant for the 
lowest excited states with S = (S = 1) for the eigen- 
values of the singlet states (triplet states) . e* , J and t^ 
are the key energy scales that determine the properties 
of the riL = 1 states. For states with purely til — 1 char- 
acter, varying e* is equivalent to varying any one of e, 
Um^ Uh, Vlm, or Vhl as given by Eq. [S] If we relax the 
approximation t^ — we need to deal with states with 
ni, ^ 1. If t^ ^ 0, then in principle e*, J and t^ are 
no longer the only energy scales - one also needs to con- 
sider the individual effects of each of the direct Coulomb 
terms. For our typical parameter set, the effect of t^ on 
the lowest excited states is quite small. This means that 
e* is still the key parameter, as it is in the t^ = case. 
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FIG. 3. Schematic of the energy and character of the rii = 1 
excited states in the regimes e* < J/4, e* ~ J/4 and 
£* > J/4, where e* = £ + Um — Uh + Vlm — Vhl- As e* 
increases, the two triplet states approach the energy of the 
singlet MLCT state. At e* = J/4 there is an avoided crossing 
and the triplet states change character, with the lowest energy 
triplet state going from predominantly LC when e* < J/4, to 
predominantly MLCT when e* > J/4. This can be seen for 
some explicit parameter values in Fig. [5] 



so long as e* does not vary too much (c/. Supp. Info. 
Figs. S7-S17).65 

Fig. [3] illustrates the character of the til = 1 singlet 
and triplet states. If (e* — J/4)/t^ is small then the 
LC and MLCT triplets will be strongly mixed (note that 
this does not require that t^ is large compared to any 
of the parameters). On the other hand, in the absence 
of interactions, the level crossing would occur at e = 0, 
the point at which the one electron H^vac) and M^vac) 
states are degenerate. Fig. 2] and Fig. [5] show that the 
n^ = 1 states are well separated from the other states 
for our typical parameter set (see Appendix for more de- 
tails). It is important to stress that a large hybridization 
of the metal d and ligand tt orbitals is neither necessary 
nor sufficient to have a triplet excited state with mixed 
MLCT and LC character. This is different to what is 
discussed in Refs. [j and |4^. It is well know n^^i^^ that 
density functional theory (DFT) tends to overestimate 
the delocalisation of electron density in excited states. 
As such the conclusions drawn from DFT with regards 
to changes in electron density, for example identifying 
transitions as MLCT or LC, may not be reliable. 

In the exact solution, neither n^ nor nu are good quan- 
tum numbers. There is not much insight to be gained 
from the analytic solutions to this Hamiltonian, so we 
proceed by using some typical parameter values (dis- 
cussed further in the Appendix) and investigating around 
these values numerically. All level crossings are avoided, 
as is clear from Figs. S] and [51 These figures also show 
that the eigenstates are almost pure basis states every- 
where except very near the avoided crossings (listed in 
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\'LC^), l^MLCT^) 
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indirect 
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indirect 


10), \'MH') 


2A = £* + 2Uh -Ul+ 3Vhm - Vhl - 3Vlm 


indirect 
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TABLE I. Level crossing points for many-body states and one electron states in the one metal one ligand model and the 
Hamiltonian element coupling them. Some of the hybridizations occur indirectly (i.e. the states are not directly connected by 
a Hamiltonian element but can be at higher orders in perturbation theory). The first block of conditions are those between 
states with a Hamiltonian element directly coupling them (coupling at first order in perturbation theory). The second block 
are those between states which have no direct coupling in the Hamiltonian (higher orders in perturbation theory). The final 
block is of the hybridization conditions for the non- interacting model {Uh = Ul ~ Um ~ Vhm = Vhl = Vlm — J — 0). 



Table H]) , since t^ and i^ are small compared to all the 
other Hamiltonian parameters. Hence, (n^) — 0,1 or 2 
for each of the singlet eigenstates, except near the avoided 
level crossings. 

Table H] lists all the (avoided) level crossings in our 
model, as well as the equivalent conditions for a non- 
interacting model. J plays a key role in determining the 
level crossing for most of the states, so the singlet and 
triplet spectra are very different. By comparison, the 
level crossing points for the non-interacting states (one- 
electron states) only depend on the one electron site en- 
ergies. Thus, singlets and triplets are degenerate in the 
non-interacting model. The differences between the mix- 
ing of the singlet states and that of the triplet states are 
apparent from Figs. mandO 

Fig. m shows for our typical set of parameters how the 
eigenvalues change as a function of A, and how the singlet 
eigenstates are dominated by one basis state except near 
avoided crossings. Note that for the typical parameter 
set no singlet avoided crossings occur for A > 2 eV, i.e. 
in the range of A values found in the ligands with appli- 
cations in OLED and OPV devices (see Appendix). Fig. 
[5] shows the character of the triplet states as a function of 
e, and how the triplet states can be strongly hybridized 
even away from avoided crossings. These results indicate 
the importance of the exchange interaction J in deter- 
mining the character and energy of the excited states. 
As is apparent from Figs. |4] and [5l for a wide range 
of reasonable parameter values the 5*1 and 5*2 states of 
our model are almost pure \-^MLCT^) and \^LC^) ba- 
sis states, respectively. As such, one can calculate the 
gap between these states to be e* + 3 J/4. With an inde- 
pendent estimate of J based on isolated ligand data (see 
Appendix A.l), one can estimate e* from the energies 
of the 5*1 and ^2 states found in absorption spectra, and 



compare it to the value found from the singlet-triplet gap 
(Eq. 17). This allows for a self-consistency check on the 
model, and an estimate of its accuracy. 



IV. EFFECT OF THE MLCT CHARACTER OF 

THE TRIPLET EXCITED STATE ON ITS 

RADIATIVE PROPERTIES 

Experimental data indicate that the emitting state is 
predominantly triplet. However, there has been consid- 
erable debate over the exact character of the state, i.e., 
whether it is LC, MLCT, or a hybrid;^^"^^ The degree of 
MLCT character in the triplet excited state has been dis- 
cussed by Yersin et al. as a key indicator of a compounds 
potential as an OLED emitter >iii^i2^ Here we examine 
how the MLCT character changes in the lowest singlet 
and triplet states, and how these changes in character 
might affect the radiative properties of the triplet states. 

In an OLED device one wants to maximize the radia- 
tive decay rate compared to the non-radiative decay rate. 
It is not possible to predict how the non-radiative decay 
rate will vary from the current model, but we can ex- 
amine how to increase the radiative rate via the increase 
of the transition dipole moment. For a triplet (phos- 
phorescent) emitter, increasing the triplets transition 
dipole moment requires significant singlet-triplet mixing 
via spin-orbit coupling, which can only occur in the pres- 
ence of the heavy transition metal core.^*' The spin-orbit 
interaction couples electrons in atomic orbitals of differ- 
ent angular momentum, allowing spin flips between these 
orbitals and hence coupling singlets and triplets. In our 
model we only have one metal orbital (which is really 
some renormalized effective orbital we have identified as 
'metal'), so we insert a singlet-triplet coupling explicitly 





e* (eV) 



FIG. 4. Singlet excitation energies as a function of the lig- 
and HOMO-LUMO gap A. The triplet excitation energies 
(dashed blue lines) are also shown for comparison. The colour 
of each curve shows how the character of the eigenstates 
changes with A. Around the example value of A = 3 eV 
there is no mixing of the lowest excited singlets, meaning 
they are nearly pure basis states. There are four pairwise 
avoided crossings and one avoided crossing involving three 
levels. These occur at the maximal hybridization conditions 
listed in Table H] Away from the avoided crossings, the eigen- 
states are predominantly the basis states for these values of 
the parameters. Note that all level crossings are avoided. This 
plot is for our example parameter values (c/. Fig. [2] and the 
appendix). Figs. S7, S9, Sll, S12, S15, S16 and S17 in the 
Supp. Info, show that for the lowest excited states, this pic- 
ture holds for a wide range of parameters around our typical 
parameter set.— 



between states with an unpaired spin in the metal orbital 
to mimic the true effect of a spin-orbit interaction. 

The triplet state \Tm) has a transition dipole moment 
Mt^ to the ground state that is given (to first order in 
the spin orbit coupling Hamiltonian Hso) by 



Mt 



E 



{Tm\Hso\S., 



Et"'' — Egn 



Msr. 



(9) 



where Et^ is the energy of the triplet state, Es^^ is the 
energy of the singlet state \Sn) with transition dipole mo- 
ment to the ground state Mgn , and the sum runs over 
all singlet states (see page 271 of Ref. [43 for more de- 
tails). This expression is only valid if 



{T^\Hso\S„) 



Ej 



'Es 



< I. 



FIG. 5. The triplet excitation energies as a function of the 
ligand HOMO-metal orbital gap e* . The singlet excitation en- 
ergies (dashed blue lines) are also shown for comparison. The 
colour of each curve shows how the character of the eigenstates 
changes with e. One can see the '^MLCT^ - '^LC^ hybridiza- 
tion point occurs at e* — J/4 = 0.25 eV as expected, while the 
^MH^ - "^LC^ hybridization occurs at very large values of e* for 
this parameter set. The location of the ^MLCT^ - ^LC^ maxi- 
mum hybridization point at 0.25 eV puts it in the likely range 
of e values, meaning that in realistic circumstances the low- 
est triplet state will probably be a hybrid state. In contrast, 
avoided crossings in the lowest excited singlet state occur at 
£ ^ 1 eV, outside the reasonable range in these complexes. 
This plot is for our example parameter values (c/. Fig. [2] and 
the appendix), varying e between and 1 eV. Figs. S8, SIO, 
S13 and S14 in the Supp. Info, show that this picture holds 
for a wide range of parameters around our typical parameter 
set, and reinforce the relevance of e* as the key parameter for 
the n_L = 1 excited states.— 



If the singlet and triplet are nearly degenerate, then one 
needs to exactly solve the entire singlet-triplet Hamilto- 
nian with explicit spin-orbit coupling. 

Fig. [5] and Figs. S2-S6 of the Supp. Info, show that 
the energy of the MLCT singlet lies between the two 
riL — 1 triplet states over a wide range of possible pa- 
rameter values, while the LC singlet is higher in energy. 
These figures also show that when the lowest triplet state 
is predominantly MLCT {i.e. e* > J/4:), and is nearly de- 
generate with the singlet MLCT state. In this regime the 
perturbation expression above is no longer valid. For rea- 
sonable parameter values the MLCT singlet is the closest 
singlet (in energy) to the lowest triplet state (see Figs. |3] 




FIG. 6. Degree of MLCT character in the lowest energy triplet 
state (TipMLCT^) as a function of e* / J. As e* / J increases, 
the MLCT state becomes favorable and begins to dominate 
the lowest triplet state. At the point e* / J = 1/4 the state is 
half MLCT, half LC [{^MLCT^\T{) = (^LC^lTi) ~ 1/^2). 
Note that around our typical choice of parameters, varying e* 
is equivalent to varying e. To give the variation in e* in this 
figure we use our typical parameter values (c/. Fig. [2] and the 
Appendix), and vary e between and 1 eV. 



and [5] and Figs. S7-S17 of the Supp. Info.).^^ If we as- 
sume that the only contribution to the transition dipole 
moment of the lowest triplet state comes from the singlet 
which is nearly pure MLCT we can rewrite Eq. IH] as 



Mti = Ml 



MLCT^ - 



(TipMLCTi) 



x{^MLCT^\Hso\^MLCT^). (10) 

Fig. inishows the amount of MLCT character in the lowest 
energy triplet state {^MLCT^\Ti) as a function of s* /J. 
We write our eigenstates as 



i^.: 



c%\0) 



c//^^^'\'MLCT') + . 



(11) 



so the i singlet states transition dipole moment to the 
ground state j^o) is 

Ms. 



E 



^4- J>' 



WW) 



(12) 



where |0),|0') are any of the basis states. We expect 
that due to the spatial separation between the metal 
and ligand orbitals, we can set (0|f|^MLCT'^) = and 
similarly for other terms involving separated excitations. 
The only terms we expect to remain are diagonal terms 



and intraligand excitations i.e. 



'LC' 



In 



the range of our typical parameter values the lowest ex- 
cited singlet state has almost pure \' M LCT^) character. 



whereas the ground state is almost pure |0), so the dom- 
inant eigenstate coefficients will be Cg and c^*^^*^^ . In 
the range of our typical parameter values, we find that 
^0 „o 



-Cg*^^*-^"^ Cg*^^*-^"^ , so we can approximate the 



lowest singlet's transition dipole moment as 



Ms. = ec^c^Jr, 

where Sr = {' M LCT'\f\' M LCT') - 
Eq. [13] into Eq. [10] we find 



(13) 



|f|0). Inserting 



M 



Ti 



,0 „0 



eCSiCSo 



{Ti\^MLCT') 



Eti — EiMLCT^ 

xSrCMLCT'lHsol^MLCT'). (14) 

ecg c'g J _^ — will vary with our parameter val- 
ues, while Sr{^MLCT'\Hso\^MLCT') is a function only 
of our basis states and the strength of the spin-orbit cou- 
pling on the metal atom. 

The radiative decay rate of the triplet, found via the 
Einstein A coefficient, is 



1 

TT 



2uj^ 



Seohc^ 



\Mj 



(15) 



where hu is the triplet-ground state energy gap and Sq 
is the permittivity of free space.— This quadratic de- 
pendence on Mti further amplifies the effects of small 
changes in e*. 

Apart from the Hamiltonian parameters, there is only 
one free parameter in the expression for the singlet life- 
time (found from Eq. [13] via the Einstein A coefficient), 
which is the 'distance' of the transition, Sr. We solve 
our model with the typical parameter values and find 
that to reproduce a singlet radiative lifetime of ^10 ns 
(see, for example, Ref. 50) we must have Sr = 20 A. 
This distance is about twice the size of a Pd(thpy)2 or 
Ir(ppy)3 molecule. In other similar charge-transfer ex- 
citations, one often finds that the geometrical distance 
between the assumed 'donor' and 'acceptor' fragments is 
not well correlated with the dipole lengthi^ In this case, 
the large value of Sr may be due to our neglect of the 
(0|f|^LC^) term in calculation of the singlets transition 
dipole moment. The approximation that there is a large 
spatial separation between the metal and ligand orbitals 
may also be the source of the discrepancy, as the effec- 
tive orbitals of our model may not be as localized as our 
labels suggest. 

We choose a reasonable value for the strength of the 
spin orbit coupling, {^MLCJ^jHsol^MLCT') - 100 
cm~^ = 0.012 eV (c/. Ref. [53) and use the value Sf 
= 20 A discussed above, and plot the calculated the low- 
est exited singlet and triplet states lifetimes as a func- 
tion of e*/J in Fig. [7] Fig. [7] shows that the triplet 
lifetime tt changes rapidly as a function of e*/J as 
long as e* /J < 1/4, i.e. the triplet state is not within 
{^MLCT'\Hso\^MLCT') of the singlet energy (the con- 
dition for perturbation theory to be valid). Changing 
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FIG. 7. Triplet radiative lifetimes tt (upper three curves, 
blue) and singlet lifetimes (lowest curve, red) as a function of 
£*/J for various values of Hso solved to first order in per- 
turbation theory. We have chosen Sf = 20 A to reproduce a 
singlet lifetime of order 10 ns. The triplet lifetime decreases 
rapidly as e* /J increases, up until the point where the lowest 
singlet and triplet are nearly degenerate (i.e. both MLCT) 
at which point the perturbative solution becomes invalid. As 
the strength of the spin-orbit coupling increases, the triplet 
lifetime rapidly decreases. We use our typical parameter val- 
ues (c/. Fig. [2] and the Appendix), varying e between and 
0.25 eV. 



Hamiltonian parameters. This implies that small changes 
to the chemistry of a complex, for example replacing a 
single hydrogen atom on a ligand molecule with a fluo- 
rine atom, could result in large changes in the molecules 
photophysical properties and, hence, the efficiency of op- 
toelectronic devices made from these molecules. Obser- 
vations of precisely this effect have been made in several 
systems, for example in a series of blue phosphorescent 
iridium complexes in Ref. Ij, and in a series of complexes 
based around N-heterocycHc cyano-substituted carbenes 
in Ref. ^Q. Ref. [2^ found that changing the ligand al- 
tered the singlet-triplet gap, and correlated this change 
in the gap with the changes in the radiative lifetime of 
the triplet state, finding the tt oc (£^7-1 — Eimlct^)'^ 
relationship we predict in the regime where the spin- 
orbit coupling can be included perturbatively, i.e., for 
e*/J < 1/4. 

In the regime e* /J > 1/4, where the lowest excited sin- 
glet and triplet are both of MLCT character, the dom- 
inant splitting between them comes from an exchange 
interaction between the metal orbital and the ligand or- 
bitals, an interaction we have neglected in this model. 
The relative size of this inter-site exchange and the spin- 
orbital coupling element will determine the excited state 
properties in this regime. Spin-orbit coupling may also 
have important effects on the non-radiative lifetime of 
the excited state, but these effects cannot be studied 
within the current model. One can expect that since 
non-radiative processes are thermally activated, small 
changes in energy barriers can cause exponentially larger 
changes in the rate. 



CHARGE INJECTION 



e*/J from 0.1 to 0.25 decreases the triplet's radiative 
lifetime by an order of magnitude. For example with 
{^MLCT^\jflso\^MLCT^) = 100 cm^^, the triplets ra- 
diative lifetime changes from 17 /iS to 1.5 fis. These val- 
ues of Tt are comparable to those found experimentally 
in various organometallic complexes, for example ~ 1 /iS 
for Ir(ppy)3 (Ref. IH), or ~ lOOus for PtOEP (octaethyl- 
porphyrin platinum(II) (Ref. |5J). 

One can understand the apparently exponential change 
in Tt as follows. We are explicitly in the regime where 
\Ti) is dominated by \^LC^). If we treat the MLCT 
component of \Ti), {Ti\^MLCT^), perturbatively we 
find it contains a factor (i^j-i — E3mlct'^)~^ = {Et^ — 
Ei-MLCT'^)~^ ■ Thus the transition dipole moment of the 
triplet contains a factor (i^yi ~ -E-ia/lctO"^! so the life- 
time varies with {Et'^- — Eimlct'^)^- Over the small 
regime we vary e*, (Et^ — Eimlct^) decreases linearly 
with £*, as seen in Fig. [5] In this small range, a large 
power will be approximately linear on a semilogarithmic 
plot. 

These results show that the excited state character 
and lifetime can be sensitively dependent on changes in 



An important problem to understand in organometal- 
lic complexes is the energetics of charge injection and ex- 
traction, and the influence of the character of the states 
involved.'^" When electrons (holes) are injected into a 
bulk sample of our organometallic complex at an anode 
(cathode), a five (three) electron state is formed on a 
complex. When these oppositely ionized complexes are 
near each other they can return to charge neutral states 
(charge recombination) with one complex in its neutral 
ground state, and the other in a neutral excited state. 
It is the photoemission from this excited state which is 
desired in an OLED device. We would like to predict the 
relative population of each emitting excited state. 

The typical lifetime of a singlet state is on the order 
of 10 ns, and for a triplet it is much longer (around 10 
/is). It is believed that in complexes such as these a 
singlet-triplet intersystem crossing can occur in 50 fsJ^ 
This is more than 5 orders of magnitude greater than 
either excited state lifetime. As such, the population of 
the excited states has plenty of time to become a thermal 
population. This means we can predict the probability of 
finding either triplet or singlet excited states in the ther- 
mal population based on the energy difference between 
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FIG. 8. Probability of injected charges forming a triplet ex- 
cited state versus e* for various values of J. This probabil- 
ity is highly sensitive to e*/J. It is clear that increasing e* 
suppresses triplet formation. The point at which the rapid 
decline from probability 1 begins is e* /J = 1/4, the point at 
which the triplet state gains large MLCT character, bringing 
its energy close to that of the lowest excited singlet state. The 
values of J correspond to the ligands ppy (J = 2.12 eV), thpy 
(J = 1.54 eV), fluorene (J = 1.17 eV) and bzq (J = 0.88 eV), 
all at 300 K. This plot was made with the typical parameter 
values (c/. Fig. [2] and the appendix), varying e from eV to 
2eV. 



those states. 

The probability of the excited state being a triplet 
state at a temperature T given by a Boltzmann distri- 
bution is 



1 + ie-^(B«-£;^)/(feBT) 



(16) 



whilst the probability of finding it in a singlet state is 
Ps = l- Pt- 

Increasing e* rapidly suppresses the triplet population, 
cf. Fig. [HI As e* increases, the probability of finding 
the excitation in the triplet state decreases, towards the 
limit of 75%. The triplet probability depends on precisely 
what ligand and metal core one uses to construct the 
complex due to the probabilities sensitivity to e* . At 
e* = 0.5 eV on the curve with J =1.17 eV the probability 
of finding the excitation in the triplet state is above 90%. 
The oft quoted proportion of 75% triplets, 25% singlets 
(see, for example, Ref. ISa ) is only reached in the limits 
£* /J -J> oo or T -J' cxo. 

In the limit t^ = we can find the lowest singlet 
(triplet) excited state by solving Eq. S] (Eq. [SJ. Thus 
we have the energy of the state with recombined charges 



E^- 



1 



3J 



- 2A + 3e* + — - Wm + Wh 



■ 6Vhl + ^Vhm 




^) +^{t"Y 



or an excited triplet 



1 



J 



E' = - 2A + 3£* - - - 2Um + 4C/h 



■Whl + '^Vhm 



Wr-i' ^'^'"^' 



The energy gap is 



E^ -E^ = - 



J+J(^*-'-i] +^i^")' 



8^ + ^] +4(t«)2 



(17) 



Thus we see again that small variations in parame- 
ter values can have large effects, in this case shifting the 
triplet probability exponentially as shown in Fig. \8\ As 
e* increases, there are two competing effects on the pho- 
toluminescent efficiency. One is the suppression of triplet 
production down to 75% seen in Fig. [H and the other is 
the rapid decrease in triplet lifetime seen in Fig. [71 The 
condition (e* ^ J) necessary for a large triplet transi- 
tion dipole moment also suppresses the probability of the 
formation of a triplet excited state. 



VI. CONCLUSIONS 

We have investigated an effective model Hamiltonian 
for organometallic complexes in electronic devices. We 
have seen that while the lowest singlet state is typically 
a nearly pure MLCT state, the lowest triplet states char- 
acter varies, changing from LC to MLCT with a highly 
hybridized region between. This variation in triplet char- 
acter is strongly dependent on the ratio e*/J. Impor- 
tantly, e* is purely a property of the complex and will 
depend sensitively on the ligand chemistry. The strong 
LC-MLCT mixing in the lowest triplet state means that 
a small shift in parameter values can cause large changes 
in the properties of the state (changing e* by a factor 
of 2 changes the triplet lifetime by almost an order of 
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magnitude). This sensitive dependence provides an ex- 
planation for the large observed changes in the photo- 
physical properties of organometallic complexes caused 
by small changes in the ligands (such as changing a sin- 
gle substituent atom on the ligand). As well as having a 
direct effect on the lifetime, the change in excited state 
energy which accompanies the change in hybridization 
causes a shift in the probability of finding the excitation 
in the triplet state. As e* increases, the triplet decay 
rate increases by orders of magnitude while the triplet 
probability decreases by at most 33%. 



ACKNOWLEDGMENTS 

We are grateful to Arthur Smith for critically read- 
ing this manuscript, and to Paul Burn, Arthur Smith, 

I 



Paul Shaw, Lawrence Lo and Seth Olsen for helpful dis- 
cussions. B. J. P. was the recipient of an Australian 
Research Council (ARC) Queen Elizabeth II Fellowship 
(project no. DP0878523). R. H. M. was the recipient of 
an ARC Australian Professorial Fellowship (project no. 
DP0877875). 



Appendix A: Estimated parameter values 

We would like to understand the values of, and rela- 
tionships between, the Hamiltonian parameters (J, A, 
t , t , Uh, Ul, Um, Vhl, Vhm, Vlm)- 

We begin with a two site extended Hubbard model, 
modeling the ligand as two sites (for example the two 
phenyl rings in a biphenyl ligand) occupied by two elec- 
trons 



H = -ti'^ (c\^C2a + cLcicr) + Ue (^i^rii; + n2^n2i) + V« ^ ni„n2a 



(Al) 



- J 



Ci^Cia- By neglecting an exchange interaction between these localized orbitals, we are making the CNDO 



with 

(complete neglect of differential overlap) approximation in this small basis seti^ By transforming to a basis of 

delocalised states H^ = ^^'^^"'^ and Lj^ = '^^■^Zf^'^ , one finds that^ 



H 



+ UninH-rnHi + nL-friLi) + Vhl {nHtnLi + nniriL-f + nH-friL^ + hhiiili) 



(A2) 



r 



where 



Uh = Ul = 



Vhl 



J^Ue- Ve, 

Ui + Ve 



Ui + Wi 



(A3) 
(A4) 

(A5) 



Ligand 


Ref. 


Si (eV) 


Ti (eV) 


J(eV) 


A (eV) 


thpy 


58 


4.08 


2.54 


1.54 


3.31 


ppy 


58 


4.99 


2.87 


2.12 


3.93 


bzq 


58 


3.57 


2.69 


0.88 


3.13 


biphenyl 


page 108 of 59 


4.33 


2.84 


1.49 


3.59 


carbazole 


page 111 of 59 


3.60 


3.05 


0.55 


3.33 


fluorene 


page 118 of 59 


4.11 


2.94 


1.17 


3.53 



1. Intraligand Exchange: J ~ 1 eV: 

The first excited singlet state energy of the two-site 
Hubbard model for the ligand is Es = Vhl + ^j and 
the first excited triplet state energy is E^ = Vhl — j- 
Thus we can easily estimate J as 



J — Es — Ej 



(A6) 



Table |TI] shows the value of J estimated in this way for 
several common ligands. This data shows that using a 
value of J ~ 1 eV is realistic for investigating our model. 
While the possible values of J vary widely, the effects of 
this variation in the lowest excited states are identical to 
variations in e* . Rusanova et al. found exchange interac- 
tions of the same magnitude (around 1 eV), having per- 
formed a semi-empirical INDO/S analysis on a series of 



TABLE II. Data for several common ligands. Experimental 
energies for the lowest visible singlet and triplet states, the 
calculated spin exchange J using Eq. IA6l and HOMO - LUMO 
gap A using Eq. IA7I 



ruthenium complexes.— They also assert that the degree 
of singlet-triplet splitting is a measure of tt backbonding. 
This feature is naturally reproduced in our model as the 
only way to have singlet-triplet splitting in our model is 
via the LC^ states, the only states directly split by J. 

Figs. S9 and SIO of the Supp. Info, show that our 
specific choice of J has no effect on the qualitative con- 
clusions drawn in this paperj^ 
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2. H-L Splitting: A ~ 3 eV: 

The energy difference between ijlij||0) and ilL!|0) 
is 2i^, and they are coupled by J/2. Thus the eigenstates 
will be separated by an energy ■\/4t^ + J^/4. Table HIl 
shows that J/2 is around 0.5 eV, while the excited sin- 
glet to ground state gap is around 4 eV. We know the 
gap \/4:tj + J^/4 will be greater than the singlet ground 
state gap of '-^ 4 eV. Since we know J/2 is an order of 
magnitude smaller than this energy, we make the approx- 
imation that the states 7J|ij| |0) and L|l!|0) are eigen- 
states (implying that A — ti). We find the ground state 
is hIhUo) with energy —A + Uh- The ground state- 
singlet gap is A5 = -^ -I- A, and the ground state-triplet 



gap is At 
find 



+ A (using Uh - Vhl = J/4). Thus we 



A 



As + At 



(A7) 



We apply this to spectral data from isolated ligands to 
estimate the value of A in Table HIl finding that A ~ 3 
eV is a realistic value to use in the investigations of the 
properties of the model Eq. [T] 

Figs. S7, S9, Sll, S12, S15, S16 and S17 of the Supp. 
Info, show that our specific choice of A within the range 
of possible values in Table HI] has no effect on the lowest 
excited states, and therefore no effect on the conclusions 
drawn herei^ 



3. Direct Coulomb Interactions on the Ligand: 

Uh ~ V^i/L i; 3 eV: 

There is an empirical relationship between the Ui on 
the localized sites and the Vi between the localized sites, 
in terms of the inter-site spacing R, 



3.5- 



I-' 



.:- 3 



2.5 




v; 



R + u; 



(A8) 



FIG. 9. Predicted Ut as a function of R with the constraint 
that we match our typical value of J = f7< — V« — 1 eV. The 
dotted vertical red lines are at the carbon-carbon bond length 
in benzene 1.4 A= 2.65 as and the distance between benzene 
ring centers in biphenyl, 3.2 A. 



We expect Uh to be in the range 3.25 eV to 2 eV (based 
on the above range for Ui). Uh will be much larger than 
0.25 eV, and hence we can approximate Uh — Vhl- This 
analysis makes it seem reasonable to choose Uh — Vhl = 
3 eV for our typical parameter set used to investigate 
the model. Ref. [l^ evaluates direct Coulomb integrals 
for a series of Ru complexes, finding values ~ 4.5 eV 
(calculated via semi-empirical INDO/S computations), 
similar to the values found from the above discussion. As 
discussed in Ref. '39', this kind of calculation is at best a 
reasonable upper bound for the value of the parameters 
in an effective low energy Hamiltonian. 

Figs. S7 and S8 of the Supp. Info, show that relax- 
ing the assumption that Uh = Vhl does not cause any 
qualitative changes in the solutions of the modeli^ 



in atomic units (page 20 of Ref. 40) . If we substitute Eq. 
^31 we find that 



Ui 



J 

2 



•II 

4 



J 

'r' 



(A9) 



Fig. [3] shows that Ui should be somewhere between 2.5 
eV and 3.75 eV, given the possible range of R and a 
typical value of J = 1 eV. It is worth noting here that a 
factor of three change in bond length is only a factor of 
1.5 in the magnitude of U . 
For J = 1 eV we have 



hence 



UH = Ui- ieV, Vhl = U^ 



Uh = Vhl + ^eV. 



reV, 



(AlO) 



(All) 



Direct Coulomb Interactions Involving the 
Metal Site: Um ^ Uh, Vhl ==; Vhm: 



Since the HOMO and LUMO are in the same location 
and have the same on-site Coulomb repulsion, we find 
that the intersite Coulomb repulsion between the ligand 
orbitals and the metal will be equal. 



V, 



HM 



V 



LAI 



Rhm 



Uh + Um 



(A12) 



We expect that R12 — Rhm (where R12 is the distance 
between the two sites of our model of the ligand, and 
Rhm is the distance between the ligand and the metal). 
As long as Um ^ Uh, we will have Vhm = Vlm — 
Vhl = Uh- This reduces the six parameters for the 
direct Coloumb integrals {Uh-, Ul, Um-, Vhl, Vhm, Vlm) 
to just two (Uh and Um)- 
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If we were to assume that Uh — Um then each four- 
electron basis state gains an energy 6Uh relative to the 
case with no Coulomb interactions. Thus it is clear that 
in this approximation the direct Coulomb interactions 
have no qualitative effect on the solutions to this Hamil- 
tonian in the four electron subspace. 



TABLE III. Values of t^ for various platinum complexes, esti- 
mated with a semi-empirical parameterization and using HO- 
MOs and LUMOs from Hiickel model calculations of isolated 
ligands using the parameterization on page 284 of Ref. |63| . 
The ligand - metal bond lengths for the complexes come from 
crystallographic data contained in the references in the table. 



For our typical parameter values, the n^ = 1 states are 
well separated from the n^ 7^ 1 states. This means that 
while we are investigating the lowest excited states (the 
iiL — 1 states) varying e* captures all the same physics as 
varying e, Um, Uh, Vlm and Vhl individually. As such, 
it is convenient to choose Um = Uh and Vhl — Vhm 
and then investigate the effects of changing e* . 



Figs. S1-S4, S6-S8, S11-S17 of the Supp. Info, show 
that varying e and Um (equivalent to varying Vlm) cause 
no qualitative changes to the solutions of the model in 
reasonable parameter ranges, ^^ We must increase Um 
more than 3 eV above the typical value before there are 
any qualitative changes to the lowest excited states which 
would alter the conclusions drawn here (see Figs. S3 and 
S4),^ 



5. Hopping Intergrals: t 



: 0.1 eV: 



Pt 


thpySa 


PPy^ 


bzqS^ 




0.09 eV 
0.11 eV 


0.08 eV 
0.11 eV 


0.08 eV 
0.06 eV 



Using the standard semi-empirical parameterization 
(for example, page 551 of Ref. |64[ ). along with 
Hiickel HOMO and LUMO orbitals of an isolated lig- 
and and experimental carbon-metal and nitrogen-metal 
bond lengths, we estimate the values of t-^ and t^ for var- 
ious ligands given in Table IIIII Note that the variation 
in the t values is almost completely due to the differ- 
ences between the Hiickel orbitals (the variation due to 
the different experimental bond lengths is ^ 1%). 



6. HOMO-metal splitting, e, is a property of the 
complex: 

e is a property only of the whole complex which is 
difficult to predict a priori. In Fig. [4] we choose a value 
of £ = 0.25 eV for the sake of concreteness. Figs. S6, S15, 
S16 and S17 in the Supp. Info, show that this choice has 
no effect on our conclusions regarding the lowest excited 
states j^ 

Figs. S1-S4, S6-S8, S11-S17 of the Supp. Info, show 
that varying e and Um and hence e* cause no qualita- 
tive differences in the solutions of the model, as discussed 
above in the section on direct Coulomb integrals involv- 
ing the metal sitei^ 
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One might be concerned that the conclusions we drawn 
from numerics at 'typical parameter values' are not ro- 
bust to variations in those parameters. Here we show 
that nothing in the lowest excited states qualitatively 
changes when we vary the parameters in a reasonable 
range. The typical parameter values around which we 
vary are J = 1 eV, A = 3 eV, e = 0.25 eV, t" = 0.1 
eV, t^ ^ 0.1 eV, Um = Ur = Ul ^ U = 3 eV, 
Vhl = Vhm ^Vlm^V^3 eV. 

Figs. S1-S4 show that the assumption Um — Ur does 
not qualitatively change the lowest excited singlet state 
until Um > 5 eV. Fig. S5 shows that the approximation 
t^ = is a valid approximation. Fig. S6 shows that 
increasing the value of e has no significant effect on the 
lowest excited singlet state. Figs. S7 and S8 show that 
the approximation V — U does not have any strong ef- 



fect on the character of the eigenstates. Figs. S9 and 
SIO show that the value of J chosen does not have any 
strong effect on the character of the singlet eigenstates, 
but does change the value of e at which the lowest triplet 
state changes character. Figs. Sll and S12 show that 
the approximation Um = Uh does not have any strong 
effect on the character of the lowest singlet state. On the 
other hand. Figs. S13 and S14 show that in the triplet 
states, increasing Um - Uh shifts the \^MLCT^)-\^LC^) 
degeneracy condition to lower values of e (equivalent to 
decreasing J).) Figs. S15 and S16 show that the exact 
value of £ chosen has no effect on our conclusions regard- 
ing the lowest excited states. Fig. S17 shows that in the 
lowest excited states varying Um is equivalent to varying 
e, as expected since e* is the controlling parameter of the 
riL states. 
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FIG. 1: This figure shows the singlet eigenstates with Um = 4 
eV, for comparison with Fig. 2, showing that there is very 
little change in the lowest singlet states. 




FIG. 2: This figure shows the singlet eigenstates with Um ~ 
4.7 eV, for comparison with Fig. 2, showing that although 
the \^LC^) and \^A4LCT^) states are mixing strongly, the 
ground state remains |0) and the first excited singlet remains 
predominantly \^MLCT^). This shows that varying Um (and 
hence e*) by 1.7 eV has no qualitative effect on the lowest two 
singlet states. 



MLCP 



W LC 



E MLCT^ 



o 



MH= 



LC= 



&: 



10.8 



I0.6 



10.4 



0.2 



12 3 4 5 6 

Eigenstates 

FIG. 3: This figure shows the singlet eigenstates with Um ~ 
5.5 eV, for comparison with Fig. 2, showing that the |0) 
and l^MLCT^) states are mixing more strongly, and that 
the \^LC^), l^MLCT^) and \^MH^) states have changed or- 
der. Even so, the lowest singlet states have not qualitatively 
changed. 
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FIG. 4: This figure shows the singlet eigenstates with Um = 
6.5 eV, for comparison with Fig. 2, showing that the |0) and 
I^MLCT^) states have changed order, with the ground state 
now predominantly \^MLCT^). This large qualitative change 
occurs 3.5 eV above our example value of \Jm, so we do not 
expect this ordering to be physically relevant. 
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FIG. 5: This figure shows the singlet eigenstates with t — 
eV (when til is a good quantum number), for comparison 
with Fig. 2, showing that dealing with blocks of constant ul 
is a valid approximation. 
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FIG. 6: This figure shows the singlet eigenstates with e = 1.5 
eV, for comparison with Fig. 2, showing where \^LC^) and 
\^ M LCT^) have changed order in energy. 




FIG. 7: Singlet excitation energies as in Fig. 3 using V = U — 
1/4 eV. The similarity between this figure and Fig. 3 shows 
that the approximation V — U does not have any strong effect 
on the character of the eigenstates. 
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FIG. 8: Triplet excitation energies as in Fig. 4 using V = U — 
1/4 eV. The similarity between this figure and Fig. 4 shows 
that the approximation V = U does not have any strong efTect 
on the character of the eigenstates. 




FIG. 9: Singlet excitation energies as in Fig. 3 using J = 2 
eV. The similarity between this figure and Fig. 3 shows that 
the exact value of J does not have any strong effect on the 
character of the singlet eigenstates, but does change the value 
of e at which the lowest triplet state changes character. 




FIG. 10: Triplet excitation energies as in Fig. 4 using J — 2 
eV. The similarity between this figure and Fig. 4 shows that 
the exact value of J does not have any strong effect on the 
character of the singlet eigenstates, but does change the value 
of e at which the lowest triplet state changes character. 




FIG. 11: Singlet excitation energies as in Fig. 3 using Um ~ 4 
eV. The similarity between this figure and Fig. 3 shows that 
the approximation Um = Uh does not have any strong eflect 
on the character of the lowest singlet state for the expected 
range of A > 3 eV (determined in the Appendix). 
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FIG. 12: Singlet excitation energies as in Fig. 3 using Um ~ 5 
eV. The similarity between this figure and Fig. 3 shows that 
the approximation Um = Uh does not have any strong eflect 
on the character of the lowest singlet state for the expected 
range of A > 3 eV (determined in the Appendix). 
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FIG. 13: Triplet excitation energies as in Fig. 4 using 
Um = 3.25 eV. This figure and Fig. S14 show that increasing 
Um - Uh shifts the |^MLCr^)-|^LC^) degeneracy condition 
to lower values of e (equivalent to decreasing J) . 
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FIG. 14: Triplet excitation energies as in Fig. 4 using Um ~ 4 
eV. This figure and Fig. S13 show that increasing Um — Uh 
shifts the \^MLCT'^)-\^LC^) degeneracy condition to lower 
values of e (equivalent to decreasing J). 
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FIG. 15: Singlet excitation energies as in Fig. 3 using e = 
eV. The similarity between this figure and Fig. 3 shows that 
the value of e has no important physical effect on the lowest 
excited singlet states for the expected range of A > 3 eV 
(determined in the Appendix). 



14 




FIG. 16: Singlet excitation energies as in Fig. 3 using e = 1.5 
eV. The similarity between this figure and Fig. 3 shows that 
the value of e has no important physical effect on the lowest 
excited singlet states for the expected range of A > 3 eV 
(determined in the Appendix). 
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FIG. 17: Singlet excitation energies as in Fig. 3 using Um ~ 4 
eV (top) and e = 1.25 eV (bottom). Both are equivalent to 
increasing e* by 1 eV. For small A values, these are very 
different. However, in the range of values expected for A 
found in the Appendix (A > 3 eV), the lowest excited singlet 
and triplet states are the same. This shows that the lowest 
excited states are only controlled by e*. If we are interested in 
higher excited states, or if we were to vary e* more, then states 
with riL ^ 1 would be important in the lowest excited states 
and varying Um would no longer be equivalent to varying e. 



